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Abstract: The dynamics of viruses are critical to our understanding of disease 
pathogenesis. Using honey bee Deformed wing virus (DWV) as a model, we conducted 
field and laboratory studies to investigate the roles of abiotic and biotic stress factors as 
well as host health conditions in dynamics of virus replication in honey bees. The results 
showed that temperature decline could lead to not only significant decrease in the rate for 
pupae to emerge as adult bees, but also an increased severity of the virus infection in 
emerged bees, partly explaining the high levels of winter losses of managed honey bees, 
Apis meUifera, around the world. By experimentally exposing adult bees with variable 
levels of parasitic mite Varroa destructor, we showed that the severity of DWV infection 
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was positively correlated with the density and time period of Varroa mite infestation, 
confirming the role of Varroa mites in virus transmission and activation in honey bees. 
Further, we showed that host conditions have a significant impact on the outcome of DWV 
infection as bees that originate fi-om strong colonies resist DWV infection and replication 
significantly better than bee originating from weak colonies. The information obtained 
fi-om this study has important implications for enhancing our understanding of 
host-pathogen interactions and can be used to develop effective disease control strategies 
for honey bees. 

Keywords: bee; viruses; Varroa; vitellogenin; temperature 



1. Introduction 

The European honey bee, Apis mellifera L. {Hymenoptera: Apidae), plays a vital role in 
agro-ecosystems by assisting in the pollination of one-third of the world's food crops. While there is a 
continuously increasing need for bees for crop pollination, the existing populations of honey bees in 
most parts of the world have undergone marked declines since the 1940s [1], threatening global 
agricultural production. An especially alarming phenomenon is the honey bee Colony Collapse 
Disorder (CCD), an event that wiped out large numbers of honey bees across the United States and 
around the world in 2006-2007 [2, 3]. A survey of microbes in CCD-affected honey bee colonies 
showed a higher incidence of viruses [4—6], raising serious concerns about risks of virus infections and 
resulting in an urgent need for effective control strategies against viral diseases to safeguard the health 
of honey bees. 

So far honey bees have been reported to be attacked by at least 19 viruses [7, 8]. Except for the 
filamentous bee virus [9], iridescent virus [10], and Chronic bee paralysis virus [11], all honey bee 
viruses reported so far contain a single stranded RNA genome encapsulated in a non-enveloped 

icosahedral capsid with a diameter of 20-30 nm. Among honey bee viruses. Deformed wing virus 
(DWV) is one of the most common and prevalent viruses infecting honey bees [3,12]. This virus, 
originally isolated from diseased adult bees in Japan [13], has a worldwide distribution and has been 
reported in Europe, North America, South America, Africa, Asia, and the Middle East [7,8,14]. DWV 
contains a positive-sense, polyadenylated, and monopartite monocistronic RNA genome encoding for a 
large polyprotein with structural capsid proteins residing at the N-terminal section of the polyprotein 
and nonstructural proteins arranged in the order of a RNA helicase, a chymotrypsin-like 3C protease 
and an RNA-dependent RNA polymerase (RdRp) residing at the C-terminal section of the polyprotein. 
Based mainly on its virion structure and genome organization, DWV has been assigned to the genus, 
Iflavirus, family Iflaviridae [15,16]. 

DWV infection is among the most common viral infections in honey bees and the virus has 
established a persistent infection in most apiaries in the world. Infection by DWV was found in 
different bee castes including queen, drones, and workers as well as in different bee developmental 
stages including eggs, larvae, pupae and adults [17]. While DWV usually persists as an asymptomatic 
infection transmitted orally between adults through trophallaxis and vertically from mother queens to 
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offspring in bee colonies [18-21], the virus can be re-activated at any time as a result of various host 
stress triggers, causing symptoms of illness to emerge in infected bees. The manifestations of the 
disease caused by DWV infection include shrunken, crumpled wings, decreased body size, 
discoloration in adult bees, and reduction in life span [13,22,23]. The severe symptoms of DWV 
infections are associated with infestations of the parasitic mite, Varroa destructor. Both laboratory and 
field studies showed that the Varroa mite is not only an effective vector of DWV by acquiring the 
virus fi-om infected bees and transmitting it horizontally to uninfected bees during feeding [24—31], but 
also an activator of the virus by causing host immunosuppression and thereby fostering viral replication 
in bees [32]. Observation of strikingly high rates of DWV infection in association with Varroa mites 
suggests that DWV is a significant factor for bee mortality and colony losses. Therefore, there is a 
critical need to understand more flilly the mechanisms responsible for inducing overt disease in 
inapparently or latently DWV-infected bees. 

Viruses have evolved persistent and acute life strategies in hosts that can change populations 
through disease-induced mortality. The dynamics of viral replication, in turn, are dependent on the 
conditions of the hosts that are regulated by biotic and abiotic factors singly or in combination. The 
levels of viral disease expression in honey bees likely result from complex molecular and cellular 
interactions between the host and the pathogen. In order to gain more insight into the incidence of 
persistent DWV infection in honey bees, we conducted a field screening to investigate the abundance 
and seasonal changes of DWV in bees under natural conditions. Using DWV as a model pathogen, we 
also conducted laboratory experiments to investigate the dynamics of the virus replication in honey 
bees following stressful manipulations. The effects of abiotic and biotic factors on the DWV 
replication were assessed by (1) experimentally incubating sealed brood at different temperatures and 
comparing the DWV levels in newly emerged bees incubated either at optimal growth temperature or 
temperature that was cooler than normal condition, and (2) experimentally challenging newly emerged 
adult bees originating fi-om strong and weak bee colonies with different levels of Varroa mites and 
comparing the DWV titer in different treatment groups. In addition, gene expression levels of 
vitellogenin (Vg), a precursor of egg yolk protein that is a plausible marker of general robustness and 
nutritional status of honey bees, were compared in bees from both strong and weak colonies to 
examine the correlation between the host health status and the virus titer in honey bees. 

2. Results and Discussion 

2.1. Persistent Infection of DWV in Honey Bees 

DWV has established a persistent infection within honeybee hosts in nature. While the vast majority 
of the DWV infected bees appeared to be asymptomatic, DWV infection was found in every apiary 
investigated and in every month of the year during the study. The percentage of pooled DWV infection 
rate was 25% and 35% for strong colonies and weak colonies, respectively, in the first month of the 
study (Figure 1). The infection rate began to rise after the spring and reached a peak of 53% and 77.5% 
for strong colonies and weak colonies, respectively, in the early winter month of December before 
declining by the end of winter. Although there was no statistically significant difference in DWV 
infection rates between the strong and weak colonies during the spring, the infection of DWV was 
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detected more frequently in weak colonies than in strong colonies after spring. The infection rate of 
DWV in bees originating from weak colonies was significantly higher than bees originating from 
strong colonies in the summer, fall, and early winter (p < 0.01). All weak colonies could not survive 
the cold weather and died before February. 



Figure 1. Seasonal activities of Deformed Wing Virus (DWV) infection in strong and 
weak honey bee colonies. DWV infection was detected in every month of the year but the 
infection rates varied from month to month during the study. All weak colonies died in 
January before the completion of the survey. 
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2.2. Elevated DWV Titers in Cold-Stressed Honey Bees 

A nest temperature of about 33-35 °C was reported to be optimal for honeybee brood development 
[33,34]. Brood incubated at 29°C failed to emerge. The reduction in temperature from 33 °C to 30 °C 
was accompanied by the addition of four extra days to the 12 day period normally required for pupae 
to emerge as adults, and a 15% decrease in total number of adult bees that emerged. Further, the 
quantification of DWV showed that the titer of the virus in newly emerged adult bees from 30 °C 
incubation was 296 ± 12.57 fold higher than bees from 33 °C incubation. The expression of 
housekeeping gene P-actin reflects the similar amount of starting material across samples 
(Figures 2A,B). 



Figure 2. The representative qRT-PCR amplification plots of Deformed Wing Virus 
(DWV) (A) and p-actin (B). 
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2.3. Elevated DWV Titers in Varroa-Challenged Bees 

Real time RT-PCR showed that DWV levels were higher in bees that originated from both sfrong 
and weak colonies after Varroa mite challenge. While DWV was also detected in some of the bees 
collected at day 0 without fransferring into rearing cups and exposure to Varroa mites, this negative 
confrol group had the lowest level of DWV infection among all examined groups and therefore was 
chosen as a calibrator. The virus level in other groups was expressed as an n-fold difference relative to 
the calibrator, the negative control group. 

For bees taken from weak colonies (Figure 3A), the DWV titer steadily increased over the whole 
experiment following the initial Varroa mite introduction and reached a maximum at day 7 
post-treatment. The highest DWV titer was observed in the group of bees with 30% Varroa mite 
parasitism at day 7 post treatment. There was a significant positive correlation between the dose of 
DWV and number of Varroa mites introduced at each time point following Varroa mite exposure 
(R^ = 0.82 at day 3 and = 0.99 at day 7; two-way ANOVA F = 5.059, p = 0.0171). 

For bees taken from sfrong colonies (Figure 3B), DWV concenfration showed an increase, and a 
positive correlation between the dose of DWV and the level of mite parasitism, similar to the profile 
observed in bees from weak colonies (R = 0.88; two-way ANOVA F = 1 1.74, p = 0.0013). However, 
following the initial increase of DWV at the early time-point post Varroa mite infestation (Day 3), the 
DWV levels declined at the later time-point of Varroa infestation (day 7). Additionally, the titers of 
DWV did not differ significantly across groups challenged by different levels of Varroa mites at day 7 
post Varroa exposure (R = 0.66), as seen in bees from weak colonies. 

2.4. Vitellogenin Expression 

The present study confirmed our previous study's findings showing that the vitellogenin transcript 
level in bees originating from strong colonies was significantly higher than bees from weak colonies 
prior to Varroa challenge [35], and further demonstrated that the gene expression of vitellogenin in 
bees was independent from the levels and time period of Varroa mite infestation and DWV infection. 
No difference in vitellogenin gene expression across groups with different levels of Varroa mites for 
bees from both strong and weak colonies at different time points post Varroa exposure was observed 

2 2 

(day 3 R =0.51; day 7 R =0.48). For the sake of clear illusfration, the relative expression of 
vitellogenin obtained at different time points post Varroa mite exposure at each category was pooled. 
While the expression of vitellogenin in each category of bees was not correlated with levels of Varroa 
mite challenges, a significant difference in the titer of the vitellogenin was observed between bees 
from sfrong and weak colonies (Figure 4) (two-way ANOVA F = 12.15; p = 0.001 1). 
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Figure 3. Relative quantitation of the DWV load in bees challenged by Varroa mites. 
(A) Bees originated from weak colonies. (B) Bee originated from sfrong colonies. Bees 
from strong and weak colonies were inoculated with variable levels of Varroa mites 
(0%, 10%, 20% and 30%) to induce the DWV infection and then compared the titers of the 
virus at different times after Varroa mite challenge. The bees collected at day 0 without 
transferring into rearing cups and exposing to Varroa mites were used as a negative 
control. While DWV was also detected in some bees in the negative control group, these 
had low levels of DWV infection. The virus level in other groups was expressed as an 
n-fold difference relative to the negative confrol group. 
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Figure 4, Relative quantitation of Vitellogenin (Vg) transcripts in bees from strong and 
weak colonies after Varroa mite infestation. 
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3. Experimental Section 

3.1. Honey Bee Colonies and Varroa Mites 

Honey bee colonies were maintained in two different USDA Bee Research Laboratory apiaries in 
Beltsville, MD. The Varroa mites were monitored by the natural drop detection method and treated as 
needed to ensure that mite levels remain below damaging levels throughout the study. The strength of 
the honey bee colonies was determined based on the number of adult bees and brood as well as the 
amount of food stored in the combs. The bee colonies that had more than ten frames covered with adult 
workers and more than six frames filled with brood and food stores were defined as strong colonies, 
while bee colonies that had a small number of foraging bees flying in and out, less than ten frames of 
adult bees, less than six combs with brood and small patches of food stores were defined as weak 
colonies (Figure 5). The colonies that were heavily infested by Varroa mites were left untreated and 
used as Varroa resource colonies. They were maintained in a separate apiary. 

Figure 5. Demonstration of comparatively strong and weak colonies in the study. 
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3.2. Field Investigation of DWV Seasonal Prevalence in Colonies 

To assess the seasonality and prevalence of viral infections of honey bees in nature, five strong 
colonies and five weak colonies were chosen for virus field screening studies. Samples of 25 bees from 
each colony were collected every month for a period of one year starting in March and finishing in 
February of the following year. All field samples were individually subject to RNA extraction and 
RT-PCR analysis for DWV and other common bee viruses (data not shown). The specificity of PCR 
products was confirmed by sequencing analysis. The monthly variation of DWV infection, in terms of 
the percentage of viral infected individuals present in each colony was recorded each month. 

3.3. Laboratory Temperature Stress Experiment 

The experiments to evaluate the virus titer change in response to cold stress were focused on bees 
originating from strong colonies since our investigation of virus seasonal activities and colony 
survivorship showed that all weak colonies could not survive the winter in Maryland. To obtain worker 
bees of uniform age, the queens from strong colonies were individually confined in a brood nest for 
24 hours to lay eggs. After brood cells were capped, the brood comb was divided into three nearly 
equal parts with each piece containing around 100 capped brood cells. The cells were transferred into 
emergence cages and kept individually in incubators set at 29 °C, 30 °C or 33 °C. The length of time 
required for pupa development and numbers of emergent adults were recorded. The emerged adults 
were removed from the cages and subjected to subsequent RNA isolation and qRT-PCR analysis for 
DWV. The procedure was repeated three times. 

3.4. Laboratory Varroa-Challenging Experiment 

The laboratory experiments were conducted during a period between spring and summer, in the 
peak brooding season when bees were actively foraging. Frames with emerging brood were removed 
individually from four strong and four weak colonies, placed in separate mesh-walled cages and 
incubated in an insect growth chamber at 33 °C, 70% RH overnight. Emerging bees were collected the 
following day and separated into two different categories based on the conditions of the colonies they 
originated fi-om, strong or weak. Each category was fiirther divided into four treatment sub-groups 
based on levels of Varroa mites introduced. Thirty bees for each group were transferred to a 
'top-feeder' rearing cup [36] and variable levels of Varroa mites were then introduced: 0% (N = 0), 
10% (3 mites per 30 bees), 20% (6 mites per 30 bees), and 30% (9 mites per 30 bees). One subgroup of 
10 newly emerged bees from strong and weak colonies were collected individually prior to exposure to 
cage isolation stress and Varroa mite challenge to serve as negative controls and defined as day 0 in 
the study. Half of the bees along with their Varroa mites were collected three days after Varroa 
exposure and the other half of the bees were collected at day 7. 

3.5. RNA Extraction 

Total RNA was isolated by placing bees in a microcentrifuge tube individually and homogenizing 
1 mL of TRIzol Reagent (RNA extraction kit, Invitrogen, Carlsbad, CA, USA) according to the 
manufacturer's instructions. Following precipitation and centrifugation, the resultant RNA pellets were 
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re-suspended in UltraPure DNase/RNase-free distilled water (Invitrogen, Carlsbad, CA, USA) in the 
presence of Ribonuclease Inhibitor (Invitrogen, Carlsbad, CA, USA). The concentration of total RNA 
was determined by measuring the absorption at 260 nm and purity of RNA was estimated by the 
absorbance ratio of 260 nm / 280 nm using spectrophotometer (Ultrospec 3300 pro, Amersham 
Biosciences). RNA samples were stored at -80 °C freezer until used. 

3. 6. Primers and Probes 

Sequences of primers and probe specific for DWV and a housekeeping gene, A. mellifera P-actin, 
were reported previously [17]. A previously published primer pair for detection of gene franscript 
levels of vitellogenin [37] was employed for measurement of general robustness of honey bees in this 
study. All the primers were sjoithesized by Invitrogen and the TaqMan probes were purchased from 
Applied Biosystems. 

3.7. TaqMan Real-Time Quantitative RT-PCR (qRT-PCR) for DWV 

The quantification of DWV titer in individual bees was performed by TaqMan real time qRT-PCR 
using a Sfratagene Mx3005P™ Multiplex Quantitative PCR System. The amplification was carried out 
using the one-step Access RT-PCR System (Promega, Madison, WI) incorporated with a 0.2 |a,M 
TaqMan probe. The amplification profile consisted of one cycle of 48 °C for 45 min followed by 
40 cycles at 95 °C for 30 s, 55 °C for 1 min, and 68 °C for 2 min. The amplification of P-actin was 
performed for all of the samples analyzed to normalize the results. The final PCR product was 
confirmed by 1.5% agarose gel electrophoresis with ethidium bromide. Negative confrols (no 
template) were included in each run of the reaction and yielded no products. 

The DWV titers in bees were determined based on the value of the cycle threshold (Ct), which 
represents the number of cycles needed to generate a fluorescent signal above a predefined threshold 
and therefore is inversely proportional to the concenfration of the initial target that has been amplified. 
The Ct values of the PCRs from each subgroup were averaged and relative quantification of DWV 
titers was performed using the comparative Ct method, also referred to as the AACj method. Briefly, 
the average amount of DWV titer of each group was normalized to internal control P-actin following 
the formula: ACt = Ctowv " Ctp.actin- The ACt value of each group was relative to a calibrator which 
represented the minimal virus level following the formula AACt = ACt(target) ~ ACt(caiibrator). The 
fold-change in the DWV concentration of each treatment group was calculated using the formula 
2-AACt comparative Ct method was validated by showing that the amplification efficiencies of 
DWV and P-actin were equal. 

The comparative Ct method was used for quantification of DWV titers. In order to qualify the use 
of the comparative Ct method for DWV quantification, the amplification efficiencies of the TaqMan 
real time qRT-PCR assay for both DWV and P-actin need to be approximately equal. The 
amplification of six five-fold dilutions of total RNA ranging from Ijxg to 0.32 ng per reaction was 
performed in triplicate and the standard curves for DWV and P-actin were generated by plotting the Ct 
value against the corresponding input RNA. A linear relationship was observed between the amount of 
input RNA and the Ct values for both DWV and P-actin. The coefficient of correlation (R ), which 
indicates the linearity of the Ct values plotted in the standard curves, was 0.9822 and 0.9572 for DWV 
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and P-actin, respectively. To control for variation in RNA samples, the concentration of DWV in all 
tested samples was normalized to the housekeeping gene, P-actin (ACj). Relative efficiencies of DWV 
and P-actin amplification were plotted as ACt = (Ct(dwv) ~ CT(p-actin)) versus the log of the 
corresponding amounts of input RNA. As demonstrated in Figure 6, the efficiency plot between log 
input RNA and ACt had a slope less than 0.1 (slope = 0.056), indicating that efficiencies of DWV and 
P-actin amplification were approximately equal (Figure 6). Therefore, the comparative Ct for relative 
quantification of DWV in bees challenged by different levels of Varroa mites was valid. 



Figure 6. Validation of amplification efficiency for DWV and P-actin by TagMan 
qRT-PCR. The amplification of six five-fold dilutions of total RNA ranging from 1 |j,g to 
0.32 ng per reaction was performed in triplicate and the standard curves for DWV and 
P-actin were generated by plotting the Ct value against the corresponding input RNA. The 
difference between the Ct value of DWV and that of P-actin (ACt) was plotted against the 
log of the corresponding amounts of input RNA. 

Validation of Amplification Efficiency for D\A/V and beta-actin 
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3.8. SYBR-Green Real-Time qRT-PCR for Vitellogenin 

The quantification of vitellogenin in bees from strong and weak colonies was performed using 
SYBR Green real-time qRT-PCR. RT-PCR reactions were carried out in a 50-\iL reaction volume 
containing 25 \iL of 2x Brilliant® SYBR® Green QRT-PCR Master Mix (Stratagene, La JoUa, CA, 
USA), 0.4 )j,M forward and reverse primers, and lug of template RNA. The thermal profile for the one 
step RT-PCR was as follows: one cycle at 50 °C for 30 minutes, one cycle at 95 °C for 10 min 
followed by 30 cycles of 95 °C for 30 s, 55 °C for 1 min, and 72 °C for 30 s. SYBR Green dye may 
cause non-specific binding because it intercalates with any double-stranded DNAs and emits 
fluorescence signals. Therefore, a melting curve analysis was performed after amplification to 
determine the specificity of the PCR products. The PCR products were incubated for 1 minute at 
95 °C, ramping down to 55 °C at a rate of 0.2 °C/s foUowed by 81 cycles of incubation where the 
temperature was increased by 0.5 °C/cycle, beginning at 55 °C and ending at 95 °C. The amplification 
of p-actin, was also performed for each samples to normalize the result. Negative controls (no reverse 
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transcriptase and no template) were included in each run of the reaction and yielded no products. The 
interpretation and validation of SYBR-Green Real-time qRT-PCR data for vitellogenin quantification 
were conducted using the comparative Ct method as described above. 

The profile of the melting curve showed the presence of a single homogeneous melt peak for all 
sample reactions, confirming specific amplification of vitellogenin primers used in the study (data not 
shown). The same method used for comparing efficiencies of DWV and P-actin amplification was also 
used for evaluation of efficiencies of vitellogenin and P-actin amplification. The relative efficiencies of 
vitellogenin and P-actin amplification by plotting ACt = Ctviteiiogenin ~ Ct p.actm versus the log of the 
corresponding amounts of input RNA had a slope less than 0.1 (slope = 0.0958) (Figure 7), indicating 
that efficiencies of Vitellogenin and P-actin amplification were approximately equal and that the 
comparative Cj method for relative quantification of vitellogenin in this study was applicable. 

Figure 7. Validation of amplification efficiency for vitellogenin and P-actin by SYBR 
Green qRT-PCR. The amplification of six five-fold dilutions of total RNA ranging from 
1 fig to 0.32 ng per reaction was performed in triplicate and the standard curves for 
vitellogenin and P-actin were generated by plotting the Cj value against the corresponding 
input RNA. The difference between the Ct value of vitellogenin and that of P-actin (ACt) 
was versus the log of the corresponding amounts of input RNA. 
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3. 9. Statistical Analysis 

All determinations were made at least in triplicate and the results, in terms of Ct values are 
expressed as mean ± standard deviation (S.D.). A two-sample Student's t-test was used to analyze for 
significant differences of DWV infection rate between bees originating from strong and weak colonies 
each month and the DWV titers between bees incubated into two different temperatures (30 °C versus 
33 °C). Differences were considered statistically significant if p value < 0.05. 

In order to obtain a distribution of the experimental data much closer to the normality [33], the titer 
of DWV as well as Vg transcripts were log transformed for statistical analysis. Normality of data was 
checked using the Shapiro- Wilk test and D'Agostino-Pearson test, while homoscedasticity was 
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checked with Bartlett's test using the statistical package Prism 5 for Mac OS X, ver. 5.0c; GraphPad 
software; San Diego, California, USA. A two-way ANOVA was used to compare genome equivalents 
of the DWV and Vg considering three factorial levels used in pairs: the Varroa mite inoculation (0, 10, 
20, 30%), the day post-treatment (3 vs. 7) and the type of colony (weak vs. strong). With the purpose to 
measure the strength of the linear relationship, the coefficient of correlation, R [38], was calculated 
between the Varroa mite inoculation and both viral copies and Vg copies. 

4. Conclusions 

DWV is one of the most common and prevalent viruses found in honey bees. Previous studies 
showed that DWV usually persists as an inapparent infection in bee colonies via both vertical and 
food-borne transmission pathways but infection of DWV can multiply to lethal levels to cause severe 
colony mortality in association with Varroa mites [39]. While significant progress has been made in 
elucidating the genomic structure and transmission process of DWV infection, very little is known 
about the mechanisms that regulate the persistent and acute infection in bees. In particular, host factors 
in the pathogenesis of DWV infection have not been extensively explored. This paper describes field 
and laboratory studies to examine the role of biotic and abiotic factors and host conditions on the 
replication and progression of DWV in honey bees. Our results clearly indicate that stress factors and 
host health conditions could have a profound influence on the djoiamics of the virus replication and the 
outcome of the host-pathogen interactions. 

The investigation of the virus's seasonal activities showed that a long-term persistent infection of 
DWV has been established in honey bees. Inapparent DWV-infected bees were readily identified in 
bee colonies examined. However, there was a variation in bees' resistance to DWV infections in 
different seasons. The infection of DWV is more common in fall and winter. This variation in 
fi-equency likely depends on the host's strength and weather conditions that modulate normal 
development and functions of the hosts. 

The observation of the significantly higher titer of DWV in bees from temperature-challenged bees 
suggests that cold-stress could not only cause a delay in normal brood development but also weaken 
the host immune response and increase susceptibility of hosts to viral infection. The weaker bee 
colonies that do not have sufficient bees to maintain a strong cluster for heat generation and retention 
create an increased vulnerability to virus infection and therefore easily succumbed to cold weather and 
collapsed rapidly in the early winter. Bee colonies with strong populations of healthy bees and 
adequate food stores have better regulation for brood temperature. However, significant temperature 
fluctuations during the winter could also alter colonies' abilities to regulate immune response and 
control infection, leading to an elevated virus titer and an increased severity of the virus infection and 
eventually impacting the bees' survival. The results from our study partly explains the high levels of 
winter losses associated with cold stress and support the previous flnding that DWV could act 
independently of Varroa mites to lead to colony losses [40]. 

It has long been recognized that Varroa mites are effective vectors and activators of several honey 
bee viruses including DWV. Not surprisingly, increasing the number of Varroa mites had the effect of 
enhancing the concentration of DWV in bees. The presence of DWV in newly emerged adult bees 
prior to the transmission assay indicated the persistence of DWV in honey bee colonies likely via virus 



Viruses 2011, 3 



2437 



vertical transmission pathway [18]. The accelerating effect of Varroa mite feeding on the rate of group 
infection and titer of DWV in individual adults observed in this study confirmed a previous finding 
that Varroa mites increase the infectivity of DWV not only by its ability to serve as a vector for 
horizontal transmission fi-om infected individuals to uninfected individuals, but also as an activator of 
the virus within the hosts. However, host conditions were also found to play a role in the outcome of 
the host-parasite-pathogen interactions. During the course of the experiment, bees from weak colonies 
were more susceptible to virus infection by having a significantly higher level of virus load compared 
to bees fi-om strong colonies. Additionally, while a continuously elevated virus titer was found in bees 
from weak colonies, there was a decline in virus titer at the later phase of the experiment in bees from 
sfrong colonies, after an initial increase in response to Varroa mite challenge. This result suggests the 
ability of host immunity to clear viruses and other immune complexes. This survival advantage 
displayed in bees originating from strong colonies indicates that host factors played an important role 
in determining the susceptibility and immune responses of honey bees to virus infection. 

There was a significant difference in the level of vitellogenin between bees from strong and weak 
colonies, although expression of vitellogenin was not reflected by any difference in the level of virus 
replication and degree of Varroa mite infestation. Vitellogenin is a yolk precursor protein and 
synthesized depending on pollen consumption in honey bees [41]. Vitellogenin is normally secreted 
into the hemolymph of honey bee queens before being transported into developing oocytes, and 
accounts for over 50% of total hemolymph proteins in queens [42]. However, vitellogenin is also 
sjaithesized by the sterile honey bee workers [43,44]. Aside from having a reproductive fimction, it has 
been suggested that vitellogenin serves important functions related to the immunity and longevity of 
honey bees [37,44^7]. ft was also shown that honey bees use vitellogenin to produce royal jelly [48]. 
Although vitellogenin is an important host component affecting robustness, it is not a component of 
the immune pathways. Therefore, the expression of vitellogenin in this study was found not to be 
influenced immediately in response to the sfresses caused by DWV and Varroa mites, ft is more likely 
that vitellogenin levels reflect the physiological status of the host and that physiological robustness in 
turn leads to a spectrum of immune responses in bees from sfrong colonies against pathogenic and 
parasitic invasion. The physiological immunity and social immunity of social insects have been 
gaining increased attention in recent years and has been suggested to be reasons for the relatively few 
innate immune genes found in honey bees compared with the fruit fly Drosophila melanogaster and 
the mosquito Anopheles gambiae [49]. In view of the data presented in this paper, it is clear that the 
physical status of honey bees has a profound effect on host immunocompetence to resist virus 
fransmission, replication, and disease progression, and that beekeeping practices that promote all 
aspects of bee health and keep queens well-mated and vigorous at all times can be a key area to target 
in protecting against pathogenic and parasitic threats. 

In sum, our study showed that honey bee host conditions and stress factors that are involved in the 
vulnerability to infection seem likely to play a key role in virus replication dynamics. The findings 
from this study have important implications in enhancing our understanding of host-pathogen 
interactions and provide an explanation for varying degrees of virulence of virus infections in the field 
and suggest that beekeeping management should focus on minimizing sfress as much as possible to 
keep colonies healthy and sfrong. Our next step will be to use genome-scale approaches to identify and 
characterize genes of the host that are actively involved in affecting the progression and outcome of 
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diseases to achieve a better understanding of the mechanisms of host responses to virus infections, 
which may ultimately lead to new tools for inhibiting virus replication through targeting the host and 
thereby mitigating viral disease impacts on honey bee populations 
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